† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant No. 11205123).
Some kinds of muscles can oscillate spontaneously, which is related to the dynamic instability of the collective motors. Based on the two-state ratchet model and with consideration of the motor stiffness, the dynamics of collective myosin II motors are studied. It is shown that when the motor stiffness is small, the velocity of the collective motors decreases monotonically with load increasing. When the motor stiffness becomes large, dynamic instability appears in the force–velocity relationship of the collective-motor transport. For a large enough motor stiffness, the zero-velocity point lies in the unstable range of the force–velocity curve, and the motor system becomes unstable before the motion is stopped, so spontaneous oscillations can be generated if the system is elastically coupled to its environment via a spring. The oscillation frequency is related to the motor stiffness, motor binding rate, spring stiffness, and the width of the ATP excitation interval. For a medium motor stiffness, the zero-velocity point lies outside the unstable range of the force–velocity curve, and the motion will be stopped before the instability occurs.
Molecular motors, or motor proteins, can convert chemical energy, typically obtained from the hydrolysis of ATP, into mechanical work and walk along filaments, which are related to many active processes in biological systems such as cargo transport, cell motility, and muscular contraction.[1–6] Some kinds of motors can walk hundreds of discrete steps along their tracks before finally dissociating, such as kinesin 1, mysion V, etc., they are called processive motors.[7–14] While other kinds of motors unbind from their tracks frequently, such as myosin II, Ncd, etc., which are called nonprocessive motors.[15–17] The nonprocessive motors cannot work alone, but evoke motion of large velocity when working with large groups, for example, the contraction of muscle is the cooperation of a large amount of myosin II motors.[15,16,18,19]
There are some interesting dynamics in the process of multi-motor cooperation.[20] The cooperation of processive motors can increase the run-length of the cargo particle significantly, but the velocity sometimes decreases with increasing motor number.[21–26] The cooperation of nonprocessive motors can generate large velocity and force (such as the macro-force of muscle is generated by a large amount of myosin II motors) and has a high mechanical efficiency.[15,16,27–30] Under certain conditions, collective myosin II motors can generate oscillations. For example, the flight muscle of some insects (e.g., bees and wasps) can generate oscillations with high frequency.[31,32] Skeletal muscle myofibrils have been shown to oscillate spontaneously, in vitro in the absence of any external driving.[33] Spontaneous oscillations have also been observed in a minimal actomyosin system, where the collective-myosin system is elastically coupled with its environment.[34] Some theoretical results show that dynamic instability would occur in the force–velocity relationship of the multiple-motor transport,[35–37] and that the motor system can oscillate spontaneously if it is coupled with a harmonic potential.[38,39]
For different kinds of muscles, their functions are quite different from each other. The flight muscle of some insects can oscillate with a frequency of about 100 Hz.[32] Whereas the main function of the skeletal muscle is contraction in one direction, and oscillations are seldom observed in the skeletal muscle. The function difference of different kinds of muscles may be related to their structure difference. In theory, the structure difference is usually reflected in some physical parameters which are important to the motion of the motor system. In this paper, we investigate the influence of some parameters on the dynamic instability of collective myosin II motors.
We base our model on the two-state ratchet model[35,40,41] and consider the motor stiffness, since the real myosin II motor has a long stalk (coiled coil) which is elastic.[42,43] The results show that the motor stiffness can affect the dynamic instability of the collective-nonprocessive-motor system. When the motor stiffness is large, dynamic instability appears in the force–velocity relationship of the collective-motor transport, and spontaneous oscillations can be generated if the system is elastically coupled to its environment.
With consideration of the motor elasticity, the two-state ratchet model is shown in Fig.
In the calculations, we adopt ΔE = 64 pN·nm, because the hydrolysis of ATP in the typical in vivo enviroment yields a ΔGATP of about 20 kBT,[46] and ΔE should be less than ΔGATP in this case (since ΔE + |ΔGexcitation| = ΔGATP in the typical in vivo enviroment[47]); and l = 10 nm, b/l = 0.2, as the muscle myosin can make a power stroke of about 5–10 nm in each dynamic cycle.[27,42,43] Here, we adopt Ω = 1000 s−1, i.e., ωoff = 1000 s−1 within the ATP excitation region and ωoff = 0 in the other regions (see Fig.
We now define the probability density pa(x, z, t) for a motor with its tail at x and head at z to be attached to the filament, i.e., the motor in state 1. Likewise, the probability density for a motor to be detached, i.e., the motor in state 2, is denoted as pd(x, z, t). The probabilities pa(x, z, t) and pd(x, z, t) obey the following coupled Fokker–Planck equations:
If the rate of linker relaxation for a detaching motor is fast compared to the rate of attachment, the head position distribution of the detached motor approximately satisfies the Boltzmann distribution[49]
The average force exerted by a single motor on the backbone is given by
We are interested in the steady-state solution of Eqs. (
In the calculation, we ignore the friction of the backbone. If the friction is considered, when the velocity is greater than zero, the externally applied force Fext would increase; when the velocity is less than zero, Fext would decrease, especially for a large |v|. But this change can be ignored if |v| is not too large, as mentioned in Section 2, the friction force is usually a few orders of magnitude smaller than Fmotor.
It is interesting that when the motor stiffness km is large, dynamic instability appears in the force–velocity relationship of the collective motor transport (in the range with ∂Fext/∂v < 0), which can be explained as follows. If a load is applied (Fext < 0), v decreases. The maximal load force that the motor system can carry is the extreme value of |Fext(v)|. If this load is exceeded, the system then changes the direction of motion discontinuously.[35] In the case of small motor stiffness, the velocity of the motor system decreases monotonously with increasing load force, and there is no instability in the force–velocity relationship, as shown in Fig.
We now define the velocity that corresponds to the extreme value of Fext(v) as vc. We find that vc can be positive (for large km, see Fig.
The ratio R almost cannot affect the sign of vc, but it can affect the value of vc (see Fig.
We now consider a collective of myosin II motors being elastically coupled with its environment as shown in Fig.
The spontaneous oscillation of the elastically coupled motor system is related to the dynamic instability of the motor system in Fig.
The oscillation frequency of the motor system is related to some parameters as shown in Fig.
Figure
Interestingly, some kinds of muscles can oscillate regularly with high frequency, such as the insect flight muscle (with a frequency of about 100 Hz). Spontaneous oscillations have also been observed in the skeletal muscle myosin system in vitro, with a frequency of about 0.5 or a few Hz.[33,34] It is believed that the spontaneous oscillations of the muscle might result from the elasticity coupling between the motor system and its environment. The titin molecules in the muscle can play the role of elastic elements.[18,47] The frequency of the oscillations can be affected by km, ωon, Δz/b, and ks, as mentioned above. Our result shows that the oscillation frequency can range from a few Hz to nearly 1000 Hz with different parameters, see Figs.
We studied the dynamics of collective myosin II motors based on the two-state ratchet model,[35,40,41] considering the motor stiffness. The results show that dynamic instability can appear in the force–velocity relationship of the collective-motor transport, and that the motor stiffness affects the dynamic instability of the collective motors.
We defined vc as the velocity which corresponds to the extreme value of externally applied force Fext(v). The sign of vc depends on the motor stiffness km. (i) For large km, vc is usually positive, and the zero-velocity point lies in the unstable range of the force–velocity curve, so the motor system becomes unstable before the motion is stopped, and spontaneous oscillations can be generated if the system is elastically coupled to its environment. (ii) For medium km, vc is usually negative, and the zero-velocity point lies outside the unstable range of the force–velocity curve, and the motion will be stopped before the instability occurs, so there is no spontaneous oscillation. (iii) When km is small, the velocity decreases monotonously with load increasing, so no instability occurs.
When vc > 0, the elastically coupled motor system oscillates spontaneously. The oscillation frequency is related to motor stiffness km, motor binding rate ωon, spring stiffness ks, and Δz/b. The frequency ranges from a few Hz (or even smaller, e.g., for smaller ks) to about 1000 Hz (or even more) with different parameters. Interestingly, the previous experiments in vitro have observed spontaneous oscillations of myosin systems with a frequency of about 0.5 or a few Hz.[33,34] It is also believed that the oscillations of the insect flight muscle (with a frequency of about 100 Hz) may be related to the spontaneous oscillations of their motor system.[40,47] Our result shows that the motor system can oscillate when km > 10 pN/nm, but the stiffness of the muscle myosin obtianed experimentally is just a few pN/nm.[27,43,51,52] It must be noted that we just investigate a motor system of a half-sarcomere and just choose a simple saw-tooth potential of W1 in our model. If one considers a more complicated and more realistic model, one may observe the high-frequency spontaneous oscillation of the muscle system even when km is just a few pN/nm. This needs to be further studied in the near future.
1 | |
2 | |
3 | |
4 | |
5 | |
6 | |
7 | |
8 | |
9 | |
10 | |
11 | |
12 | |
13 | |
14 | |
15 | |
16 | |
17 | |
18 | |
19 | |
20 | |
21 | |
22 | |
23 | |
24 | |
25 | |
26 | |
27 | |
28 | |
29 | |
30 | |
31 | |
32 | |
33 | |
34 | |
35 | |
36 | |
37 | |
38 | |
39 | |
40 | |
41 | |
42 | |
43 | |
44 | |
45 | |
46 | |
47 | |
48 | |
49 | |
50 | |
51 | |
52 | |
53 |